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ABSTRACT 

In the course of our VLTI young stellar object PIONIER imaging program, we have identified a strong visibility 
chromatic dependency that appeared in certain sources. This effect, rising value of visibilities with decreasing 
wavelengths over one base, is also present in previous published and archival AMBER data. For Herbig AeBe 
stars, the H band is generally located at the transition between the star and the disk predominance in flux for 
Herbig AeBe stars. We believe that this phenomenon is responsible for the visibility rise effect. We present a 
method to correct the visibilities from this effect in order to allow "gray" image reconstruction software, like 
Mira, to be used. In parallel we probe the interest of carrying an image reconstruction in each spectral channel 
and then combine them to obtain the final broadband one. As an illustration we apply these imaging methods 
to MWC158, a (possibly Herbig) B[e] star intensively observed with PIONIER. Finally, we compare our result 
with a parametric model fitted onto the data. 
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1. INTRODUCTION 

The processes that lead to the formation of exoplanets are important to understand. Stars form after a collapse 
of a giant cloud of dust and gas. After a million year, a protoplanetary disk is forming around the star, believed 
to be the birthplace of planets. 

A young star is surrounded by an active environment with which it interacts. Accretion disksP inner gaseous 
disks, 2 4 infalling envelop renmands, wuukP^land jetdHE3 are the main components of such environments. The 
complexity of physical phenomenon at play requires direct observation at the astronomical unit (A.U.) scale. 
Optical interferometry is able to bring such informations, because it can observe both in the near infrared, 
where the hot dust and hot gas nearby the star are emitting, and resolve the first A.U., which correspond to 
milliarcsecond scale at the distance of star formation regions. 

Interferometry consists in combining the light of 2 or more telescopes in order to measure the complex degree 
of coherence. For that purpose, the interferometer measures interference fringes. The amplitude of the fringes 
yields the norm, and its position the phase of a complex quantity called visibility V(u, v). Thanks to the van 
Cittert-Zernicke theorem we know that the Fourier transform of the visibilities in the Fourier Plan (u, v) gives us 
the intensity distribution I(x, y) of the source. Unfortunately, in the near infrared (NIR) the atmosphere blurs 
the phases of the visibilities. In practice, there are two interferometric measurements : the squared amplitude 
of the visibilities V 2 is the first one. The second measurement of optical interferometry is the bispectra^ and 
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Figure 1. Sketch representing three Auxiliary telescopes from the VLT Interferometer. Each telescope have a phase 
coming from the atmospheric turbulence. This phase will blur the astrophysical phase. In order to retrieve the physical 
information, we will compute the closure phase (Eq. [4] and Eq. |S|. 

particularly their phase. A bispectrum is a product of visibilities on three different baselines forming a triangle. 
If we consider three telescopes A, B and C (see Figjl]) the bispectrum will be : 

bispectrum : VabVbcVca (1) 

The phase of the bispectrum is : 

VAaVbcVcA = |^AB|e 0AB |y B c|e* BC |y C A|e* CA (2) 

4>bispectrum = 4>AB + 0BC + ^CA (3) 

Since we do not have the phases because of the atmosphere turbulence, each telescope have residual phase. 
However, the bispectrum phase is astrophysical because this quantity cancels the turbulence effects : 

J. /data i j.turb J,turb , idata , J,turb J,turb , /data . J,turb J> tur t> I A\ 

(Pbispectrum — 4>AB + <P B ~ <P A + VBC + PC ~ <Pb + PCA + 9 A ~ PC l 4 ) 

and finally : 

^bispectrum — <I>A%? + ^B^c" + ^CA" (5) 

For each pair of telescopes, one can make multiple measurements (in the (it, v) plane) thanks to the number 
of spectral channels (since the spatial frequency (/ = \fu 2 + v 2 ) is equal to the baseline length (B) divided by 
the wavelength (A) : / = B/X). This is called the spectral super synthesis. We noticed in several datasets that 
the visibility is higher at short wavelengths. If we plot the squared visibilities V 2 in function of B/X, we can 
see (Fig. [2| that the rising curve of visibilities per base is not fitting the general trend of the data for different 
baselines. First, it was seen in AMBERpS] data, but it was considered as an instrumental defect. Now, the same 
effect has been observed with PIONIER 1 ^ . We try to explain this effect astrophysically, claiming that the image 
of the object is varying through the different spectral channels inside the same spectral band, and we propose 
three techniques in order to take it into account and to be able to reconstruct images. 

These methods will be applied to an astrophysical object. They are useful to analyze MWC158 (also known 
as HD50138). This star is a Be star known to have the B[e] phenomenon and presents a strong variabilityS^HIZl 
which complexify the evolutionary stage identification of the source. Its distance is poorly constrain (d = 
500pc± 150pcP). 

In section [2] we will describe the chromatic effects in the visibilities and the section |3] will show the different 
methods to deal with them. Finally we will apply them to the astrophysical case of MWC 158 in the section [3] 
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Center : PIONIER 13 squared visiblities. Right : 
For PIONIER data the color is in function of wavelength (blue : short wavelengths, red : long 



2. CHROMATISM 

Since interferometric instruments with spectral dispersion exist, we need to take into account the flux variations 
with the wavelength in order to correctly analyze the data and have access to the spectral super synthesis. In 
the case of Young Stellar Objects (YSOs), we noticed that the visibilities have a strong spectral dependence such 
as the geometrical shape of the object could not explain it. For Herbig AeBe star, the chromatic effect explained 
in the section [l] exists typically for the Near Infrared interferometry. In the following, we explore the possibility 
that this effect is caused by a different spectral index between the central star and its surrounding media. 

2.1 Modeling the effect 

In order to confirm that, we made a simple model with a central star and its dusty disk. 

2.1.1 The star 

In our model, the star is considered to be unresolved. This hypothesis is justified for the young objects we are 
looking at. If we suppose a young Herbig star at a distance of more than 200pc and with a radius of 5 solar 
radii, then its angular radii will be O.lmas. For 100m baseline, its visibility will be V = 0.9986. To simplify our 
model, we assume V s t a r = 1< 

For the star we have 3 parameters : the radius (i?*), the distance (d) and the temperature (T*). If we asumme 
a Herbig AeBe star with a temperature of 12000K radiating as a black body, we know that in NIR we will look 
on the Rayleigh- Jeans regime of a black body (see Fig. [3]). That means that the spectral curve is proportional 
to a power law : F^ tar cx A~ 4 . 

2.1.2 The disk 

The disk model is simple : it is a geometrically thin optically thick passive disk. Its temperature is a function of 
the radius : 

T(r)=T(r Q )(^-y q (6) 

with : 

9=1 (7) 



see references 20} 21 



The disk will be sampled on several rings, each ring having its own temperature as a function of its distance 
to the star. The other geometrical parameters are the inclination (i), the inner and outer rims radii (Rim Rout), 




Figure 3. On the left, we show the location of the PIONIER spectral channels on the SED of the model. We can see 
that they are located at the crossing between the stellar and the dust fluxes. On the right, we can see that the chromatic 
phenomena is reproduced. 



and the temperature of the inner rim (T in ). The flux of each ring will be a black body at the temperature of the 
disk. The ring visibility is defined as follow^! ; 



V rmg = J Q (2*r-) (8) 



To obtain the Fourier Transform of the disk we have to add the flux of each ring and sum every contributions : 

V dlsk = Jo(27rr i -)2 7 rr iJ B A (Ti)dr i (9) 

In the results shown in Fig. [3] the chromatic effect which tends to look like the data shown in Fig. [2] The 
visibilities have the same behavior than the data. We can conclude that the chromatic effect is not instrumental 
but astrophysical. 

The effect is dominated by the flux ratio which is changing through the different spectral channels. If we 
compute the total correlated flux we have : 

V tot (B/X)F tot (X) = F.{\) + V disk (B/X)F disk (X) (10) 

with : 

F to t(A) =F,(A) + F disk (A) (11) 

If we introduce the stellar to total flux ratio /*, we obtain the mathematical description of the chromatic 
phenomena : 

V tot (B/X) = MX) + V disk (B/X)(l - /.(A)) (12) 

with 

AW " I® (13) 

In the next section we will discuss the different methods to overcome the chromatic effect. 



3. METHODS 



Our goal is to be able to analyze chromatic data. We developed three complementary methods to do that : 
gray image reconstructions, data modification and parametric fit. The first two methods are based on image 
reconstruction and the last one is model fitting. We are mostly interested in the disk around the star and we are 
looking for informations on the resolved geometry and the strength of the chromatic effects. 

3.1 Image reconstruction per spectral channel 

Once we are aware of the chromatic effect, one can make image reconstructions selecting only one wavelength 
per reconstruction (see Fig. EJ). In that case the gray image reconstruction is justified. The technique is to have 
one image per wavelength and to stack all the images in order to have the final broadband one. 




Figure 4. An image reconstruction is made for each spectral channel of the instrument. Then all the images are stacked 
together in order to obtain the final image. 

The presence of various components of different spectral indexes prevents from using a gray emission approx- 
imation in the image reconstruction process. As a consequence, since we need to work on a per-spectral-channel 
basis the (u,v) coverage quality is severely affected. 

3.2 Modification of the data 



We want to have access to the disk visibilities. From the Eq. (12), if we know the SED and then the stellar flux 



ratio /*(A) and its variation through the wavelengths, we can compute the disk visibilities as : 

v^ Bm = v ^m=m (14) 



We can apply the modification described in Eq.(14| to one of the interferometric measurements which is the 



power spectrum (VV* = \V\ ). In summary, our measurements are |Vd a ta| an d we w & n t to recover \Vd 



Using the Eq( 14 ) , we have : 

VUBIX) = ( N^^W- '•(*> )■ (i5) 

One of the problem is the value that we take for \/|Vdata(-S/A)| 2 ; we must choose between the positive (phase 
</> = 0) and the negative one (0 = it). But it could be solved analyzing more precisely the data and other 
interferometric observables like the phase of the bispectrum (also called the closure phase) . 

It is not possible to retrieve the bispectrum phase of the dust from the data because we are loosing the phase 
of each pair of telescopes (see Eq. [I| . The equations lead to a solution where we need the phase^ . 



3.3 Parametric model 



In this section we have attempted to model the object. The model is geometrical and includes the chromatic 
effect as described in the sections [T] and [2] Our model is composed of multiple components and was developed 
when chromatic data was fitted. 

3.3.1 Geometric part of the fit 

The first component of the model is an unresolved star (a dirac in the image space) which can be shifted compared 
to the image photo center (that will produce a rise of closure phases). The second component is a ring. In the 
Fourier space the ring is defined as in Eq. (JsJ) but using \/u 2 + v 2 for the spatial frequencies (B/X) and their 
orientations that we want to solve. 

In order to be able to have a Position Angle (P.A.), which is defined from the North to the East, and a 
inclination (i) we will modify the wv-plan for the extended component as follows : 

Wring = u cos P. A. + v sin P. A. (16) 

fring = (—u sin P. A. + v cos P. A.) cos i (17) 

One of the parameters of this shape is the ring radius r. But this will define a ring with a infinitely small width. 
In order to have a Gaussian width we have to convolve the ring formulae by a Gaussian, in other words, to 
multiply the visibility of the ring by the visibility of the Gaussian function with the correspondent width w : 



-(™f) 2 

^gaussianring — King exp — — — (1^) 



Once we have the Gaussian ring, we will add some azimuthal modulations of the ring intensity to be closer to 
the physics of an inner rim. The modulations are functions in cosinus and sinus of the azimuthal angle (a) of 
the ring which starts at its major axis. We have included two sorts of modulation : one on 2 tt (ci, si) and the 
second on tt (c2,S2). These modulations are defined as : 

-^tot = -fring(ci cos a + Si sin a + C2 cos 2a + S2 sin 2a) (19) 

In the visibilities that gives : 

P> B 

Vtot — King — iJi{2ir— r)(a cosa + s\ sin a) — J 2 (27r— r)(c2 cos 2a + s 2 sin 2a) (20) 
A A 



The set of parameters defines geometrically a Gaussian ring that can be fitted to the data. 

To complete the model, we can add a second Gaussian ring or a Gaussian function. The Gaussian function 
is defined with its own position angle (P.A. ganss ) and inclination (z gaU ss)- Its visibility is defined as follows : 

V gauss = e"" "™ T> (21) 

with r gauss being the Half Width at Half Maximum (HWHM) of the Gaussian function. 

We have a model with three components : the star, the Gaussian ring and a second Gaussian ring or a 
Gaussian function. 

The total visibilities are depending on all these components weighted by their flux. 



3.3.2 Modeling the chromatism 

To obtain the model visibilities we use the linearity property of the Fourier transform. 



^totHot = F*V* + FiVi + F 2 V 2 (22) 



The fluxes are the ones recieved by the interfcromctric instrument. So we have : 



F= I F x d\= / F v dv (23) 



We will use the approximation that the channel spectral width is constant and that the flux is constant in one 
spectral band. The flux is then equal to the value of F\ at the central wavelength of a spectral channel. The exact 
value of all the terms are described in the Appendix [Aj From the Eq. ( 22 1 , we see that we can determine a flux 



ratios at one wavelength and to deduce the ratios on the other wavelength by the laws that we assume for each 
component. PIONIER is operating in the NIR in the H band. At this wavelength, we can assume that Herbig 
stars are in their Rayleigh- Jeans regime. That means that their flux (F\)is proportional to the wavelength at the 
power of —4. The laws for the environment are more difficult to find. We can fit a power-law in wavelength or 
to a black body variation if we are resolving a thermal emitting region. Since the dust temperature is supposed 
to be below 2000KEED , we can assume that it is in its Wien regime. Then if we assume black body regimes we 
obtain : 

V^BA)^f^h- 4 + f 1 §^VABA) + f 2 §^V 2 (B,X) (24) 
with / the flux ratios at Aq, T the temperature of a component and 



2hc\- 5 , , 

B{\T) = rc (25) 

is the black body function, with h the Planck constant, c the light speed, and kg the Boltzmann constant. 

The variations of the flux ratios through the observational band will build the chromatic effect that we want 
to take into account in our fit. 

Once we get all our tools to investigate data with chromatic effect, let us apply them on an astrophysical 
case : MWC158. 



4. THE CASE OF MWC 158 

The interest on this object came with the data we get with PIONIERS1 a 4 telescopes interferometric, visitor 
instrument operating at the VLTI and which observe in the H band. 

4.1 Image reconstructions 

We were interested into this data (see Fig. [2]) because it shows clearly signs of chromatism. As the (u, i>)-plan is 
sufficiently covered we can reconstruct images. We use the Mira algorithm^ , but as many image reconstruction 
algorithms it does not take into account the chromatism. Since it extrapolate the Fourier space, the chromatism 
makes him extrapolate badly and many artifacts appear. We then use the monochromatic reconstructions per 
spectral channel. We also use the visibility correction, not modifying the closure phases (they should be stronger) . 
The results are showed on the Fig. [5j 

We can see that there is a second resolved component. We can also see the orientation of the smallest 
extended component. Both of the reconstruction methods shows similar patterns. That brings us to the idea to 
fit two extended components. 



Sum of images through spectral channels 




RA (mas) relative a (milliarcseconds) 

Figure 5. Left : the stack of image reconstructions per spectral channels. Right : image reconstruction after modifying 
the squared visibilities. 

4.2 Parametric fit 

The fit bring us an idea on the geometry and the light emission from the source but with a strong a priori which 
is the model we want to fit. That is why we took the geometries suggested by the image reconstructions. We 
can see a central extended part which is composed of the star and its environment which seems to have a P.A. 
and an inclination. 

We have begun the fit with one extended component which is a Gaussian function or a Gaussian ring. 
Both of the fits gave us the more or less the same inclinations and P.A. which are consistent with the image 
reconstructions. But the data was not entirely fitted : the short baselines indicates that there is a more extended 
component as showed by image reconstruction. We then add another component to our fit. In order to fit the 
strong closure phase signal we add azimuthal modulation to the ring. It appears not to be sufficient, and the 
best fit was to shift the central star. It is the only solution to fit the closure phases. 

In the end, and adding the different parameters, we ended with 15 parameters and a x 2 of 3.5. In the current 
state of the data processing and interpretation, we believe that the best fit is presented Fig. [6j The parameters 
are on the Table [I] We can see that the best fit is done with two Gaussian rings. 

The geometrical fit suggests a star, with a relatively close Gaussian ring (radius of 1.5 mas) with a lot of flux 
(py 60%). We interpret that as the resolution of the inner rim of the dusty disk. Its azimtuhal modulation is 
strong in the semi minor axis direction which leads us to deduce that it is due to the inclination. Moreover, the 
star is shifted towards the most brilliant part of the inner rim. It indicates that the inner rim has a non-negligible 
height. The outer ring suggests the continuation of the disk, or a part of the disk which is not self shadowed, or 
a halo. The constrains are poor so we can not conclude on its origins. 

The results are shown in Table [2j They are closed to the images get by reconstruction. The results are also 
consistent with that found with previous observations^!^ . The authors found similar P.A. with close values of 
the inclination of the most luminous extended object. 

The fit of the chromatism, indicates us a black body temperature of the inner rim of Ri 1500K (see Table [TJ. 
This is approximately the dust sublimation temperature found in the litteraturdHEH . 

5. CONCLUSION AND PERSPECTIVE 



The chromatic effect due to the flux predominances of two objects of different sizes is well understood and can be 
used in order to find astrophysical information of the object. In the case of Herbig AeBe stars we are able to have 




Figure 6. The best fit results are presented. Left : The image corresponding to the best fit. Center : The fit on the V2. 
Left the fit on the Closure Phases. 



an approximation of the temperature of the components. If the chromatic information is given, we can perform 
gray disk image reconstructions. They are contributing to the astrophysical analyze of the object because they 
shows the P.A. and the inclinations of the disk. Moreover, in the case of MWC158 it brought us the idea of the 
second extended component. By the fit we were able to find the inner rim radius and its temperature and to 
compare what we found with the data from photometry. We bring the first confirmation of the dust sublimation 
temperature at the inner rim. The information taken from the NIR interferometry and the chromatic effect 
argue in favor of a young nature of MWC 158. 

The main challenge is to be able to make chromatic Young Stellar Objects image reconstructions keeping 
the super spectral synthesis and without information on the total flux variation. One of the thing which is in 
process of testing, is the adaptation of the Mira algorithm to the case of young stellar object. The "gray" mira 
free parameters are the image pixels intensities. If we define the image as the image of the dust at Ao, the start 
can be represented by a dirac at the center of the image. Hence, we can put the stellar flux and a stellar relative 
spectral power law as additional parameters to the fit. Since the regularization will tend to smooth the Fourier 
plan, the algorithm will favors the added parameters to fit the fixture. The first tests was realized on MWC 158 



Table 1. The parametric fit results. The acronyms are described in the section [373] 
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Table 2. The previous interferometry results on MWC 158. Some results were complete on instruments watching at longer 
wavelengths (10.7 or 2.2 pm). The P. A. are consistent and the inclinations i also. 
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Figure 7. Image corresponding to a polychromatic reconstruction with the mira algorithm. 



(Fig. [7]). The result found a star with a flux ratio at Ao of 10.7% and a relative spectral slope of = —4.93, 
which corresponds to a black body temperature the disk image at 1.65/im of 900K. 



APPENDIX A. ANALYTIC FLUX FORMULAE 

In this appendix we will present the exact equations for the flux received by each spectral channel in order to be 
able to compute the good model of flux both for the disk and the star. 

The flux received by the instrument is : 



F= / F x dX 



(26) 



with F is the integrated flux received by the instrument, A the wavelength. 
A.l The star 

For the star we are in the Rayleigh- Jeans side of the black body flux curve, so we have : 



A+^ 



2k B cTX~ dX 



A-%^ 



(27) 



with ks the Boltzmann constant, c the light speed, T the temperature of the star. 

Since we are interested only in the evolution of the fluxes in through the spectral channels we need only : 



K, oc 



X dX = 



A^ 
3 



(28) 



A. 2 The disk 

For disk will be in the Wien regime of the blackbody curve. So we will have : 



Fdisk 



he 

2hc 2 X- 5 exp(-—-) dX 

AA KgAl 



(29) 



with h the Planck constant. 

Without the constant parts it would be : 



r-A+%^ 



A-°exp(--M dX 

AA A 



(30) 



with C\ — y^t- In order to have the analytical solution of it, we have to integrate it by parts. The first step 



is : 



Fdisk oc — ( 



A+- 



3A- 4 cxp(-^) dX + A) 



(31) 



with A = [A- 3 cxp(-^)]^A. 
Then : 



/ 3A- 4 cxp(-^) dX= ^-(2 
A 



A+4^ 



A- 3 cxp( — -i) dX + B) 

AA A 



(32) 



with B = [A- 2 exp(-^)] 



A+^ 



If wc compute 



with C = [A- 1 ea;p(-^)]^ 
Finally : 



AA 



A 2 cxp(-^) <£X = 

A _AA A Gl 



exp(-^-) 



(34) 



With all these computations the final result for the disk flux is : 



2hc2 I A 3 ( n 2 

-c^[ A+ cX B+ c[ 



cx p(-^) 



A + 



_\ A 



(35) 
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